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This study is devoted to a detailed numerical testing of the analytical results obtained recently 
for the M t -scaling of the parameters of the Bose- Einstein correlation functions. Numerical testing 
of analytical results for the spectrum is also performed, since the geometrical sizes are folded into 
^H. the momentum distribution in the scaling limiting case. The analytical results are shown to be valid 

pH | in a wider rapidity-range than thought before. 



1. Introduction 

Recently there has been much interest in the measurement and the calculation 
of the parameters of Bose-Einstein correlation functions (BECF-s) and those of the 
invariant momentum distributions (IMD-s) for rapidly expanding systems with flow 
and temperature profiles. Such systems are expected to be formed in high energy 
heavy ion reactions. The quality and the amount of available data improved drasti- 
cally mainly due to the efforts of NA35, NA44, NA49 and WA93 collaborations at 
CERN H, |, |, f§. Data from the dedicated HBT experiment NA44 indicated an 
unexpected, scaling behavior for the parameters of the BECF-s |§ . Within the errors 
of this measurement, the R S ide, Rout and R[ ong parameters of the BECF for S + Pb 
200 AGeV central reactions turned out to be equal and all were found to scale simul- 
taneously, proportionally to 1/y/Mf, where M t is the transverse mass of the identical 
particle pair. This scaling relation is found to be independent of the particle type at 
the present level of experimental precision. 
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This new generation of HBT data triggered a burst of activity mainly from the 
Budapest §, §, g g g [TOf, Kiev Q 0, |IJ, Marburg 0, |15[ and Regensburg g 0, 
0, |T^j groups, aiming at interpreting the data in terms of collective, hydrodynamic 
behavior. The results from Los Alamos group pOf also indicate collective behavior 



caused by re-scattering of secondary particles in the framework called RQMD event 
generator. These results are also supported by simulations of re-scattering effects 
in a hot expanding gas of hadronic resonances pi . It became clear, that Bose- 
Einstein correlations are in general not measuring the whole geometrical sizes of big 
and expanding finite systems, neither in the longitudinal [22] nor in the transverse 
and temporal directions, since the expansion may result in strong correlations between 
space-time and momentum space variables not only in the longitudinal, but in the 
transverse and temporal directions, too, see ref. || and references therein for a more 
detailed account on this topic. 

Where have all the geometrical sizes gone? One can show[|, H|, that they are 
disguised in the invariant momentum distribution of the bosons in case they cancel 
from the radius parameters of the Bose-Einstein correlation function (BECF). 

We briefly review here the analytical results presented in refs. |5|, |6], 0, |9], |K| which 
are tested in a detailed numerical analysis in the subsequent parts. The review of 
the analytical results mainly follows the lines of ref. [IIJ, but the present review is 
more technical than that, since we discuss here also the details of the model emission 
function and that of the IMD. Ref. [|l(| is recommended for the illustration of the 
physical ideas, which determine the M t dependencies of the parameters of BECF-s. 



2. Wigner Function Formalism 

The two-particle inclusive correlation function is defined and approximately ex- 
pressed in the Wigner function formalism as 

r-IAL. ^ (n(n-l)) JV 2 ( Pl ,p 2 ) \S(Ak,K)\* 
C(Ak-K) = w K( - m - ) ^ l+ § ID 



In the above line, the Wigner-function formalism [BaL |23] , Pq[ is utilized assuming fully 
chaotic (thermalized) particle emission. The covariant Wigner-transform of the source 
density matrix, S(x,p), is a quantum-mechanical analogue of the classical probability 
that a boson is produced at a given x^ = (t, r) = (t,r x ,r y ,r z ) with p M = (E, p) = 
(E,p x ,p y ,p z ). The auxiliary quantity S(Ak,K) = J d A xS(x, K) exp(iAk ■ x) appears 
in the definition of the BECF, with Ak = P\ — Pi and K = (pi + p 2 )/2. The single- 
and two-particle inclusive momentum distributions (IMD-s) are given by 

, N E da , E\Eo da , 

Ni(p) = — — = S(Ak = 0,p), and JV 2 (pi,p 2 = , , , (2) 

o'totdp a tot dpirfp 2 
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where a tot is the total inelastic cross-section. Note that in this work we utilize the 
following normalization of the emission function^: J —= r d A xS(x 1 p) = (n). 



3. Effects from Large Halo of Long-Lived Resonances 

If the bosons originate from a core which is surrounded by a halo of long-lived 
resonances, the IMD and the BECF can be calculated in a straightforward manner. 
The detailed description is given in ref . [[J] , here we review only the basic idea. 

If the emission function can be approximately divided into two parts, representing 
the core and the halo, S(x; K) = S c (x; K) + Sh{x\ K) and if the halo is characterized 
by large length-scales so that Sh(Q m i n ', K) « S c (Q min ; K) at a finite experimental 
resolution of Q m in > 10 MeV, then the IMD and the BECF reads as 

iVx(p) = iV liC (p) + iV(p), (3) 
I SJAk,K) I 2 

C(Ak;K) = 1 + A ' c } t ' /J , (4) 



where iVi^p) stands for the IMD of the halo or core for i — h,c and 



A* = A* (if = p) 



iVi(p) 



(5) 



Thus within the core/halo picture the phenomenological A* parameter can be obtained 
in a natural manner at a given finite resolution of the relative momentum. This 
parameter has been introduced to the literature by Deutschmann long time ago|26| . 



See ref. ]7| and references therein for a more detailed account on the origin of this 
parameter A*. In the core/halo picture, the effective or measured intercept parameter 
A(p) can be interpreted as the momentum dependent square of the ratio of the IMD 
of the core to the IMD of all particles emitted. 



4. General Considerations and Results 

We are considering jets in elementary particle reactions or high energy heavy ion 
reactions, which correspond to systems undergoing an approximately boost-invariant 
longitudinal expansion. For expansions fully invariant for longitudinal boosts, the 
emission function may depend only on such variables, which are invariant for lon- 
gitudinal boosts. These are defined as r = Jt 2 — rj = 0.51n[(£ + z) / (t — z)], 

m t = yjE 2 — pi, y = 0.5 \n[(E+p z ) / (E—p z )\ and r t = ^jr 2 x + r|. For finite systems, 
the emission function may depend on rj — too, where y$ stands for the mid-rapidity. 
Approximate boost-invariance is recovered in the | t] — y |<< Ay region, where the 
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Figure 1: Emission of particles with a given momentum is centered around r s and r) s 
for systems undergoing boost-invariant longitudinal expansion, as indicated by the 
shaded area. 



width of the rapidity distribution is denoted by Ay. In terms of these variables the 
emission function can be rewritten as 

S c (x;K)d 4 x = Sc,*(T,r),r x ,r y )dTTodr]dr x dry. (6) 

The subscript * indicates that the functional form of the source function is changed, 
and it stands for a dependence on K and y also. 

In the standard HBT coordinate system j27| , the mean and the relative momenta 
are K = (K , K out ,0, K L ) and Ak = (Qo, Q out, Q side, Ql)- Note that the side compo- 
nent of the mean momentum vanishes by definition^, [|. Since the particles are on 
mass-shell, we have = K-Ak = K Q — KLQL — K out Q out . Introducing (3l = Kl/Kq 
and p ou t = K out /K , the energy difference Qo can thus be expressed as 

Qo — PlQl + PoutQout- (7) 

If the emission function has such a structure that it is concentrated in a narrow 
region around (t s ,7) s ) i n the (r, rj) plane, then one can evaluate the BECF in terms 
of variables r and rj by utilizing the expansion 

Ak-X = Q t- QoutTx ~ QsideTy - Q lt z - (8) 
QtT Qout^x Q sidefy -Q v r s {v-Vs)- (9) 

The coefficients of the r and the t s (i] — rj s ) are new variables given by 

Q T = Qo cosh[i] s ] - Q L sinh[?7 s ] = ((3 t Q out + (3 L Q L ) coshfr^] - Q L smh[r) s ], (10) 
Q v = Ql coshfr?,,] - Q sinhfr?,,] = Q L cosh[r? s ] - ((3 t Q out + PlQl) sinh^]. (11) 

In terms of these new variables the BECF reads as 



S(Ak, | r | S C ^(Q T , Qr\, Qout, Q 



C(Ak- K) ~ 1 + 1 ~y-'"' ' ~ i + XJK) 1 ~^^y^",^,~*^e) , (12) 
l^(0,^)| 2 |5 C ,,(0,0,0,0)| 2 1 ; 

At this level, the shape of the BECF can be rather complicated, it may have non- 
Gaussian, non-factorizable structure. Gaussian approximation to eq. ( TT^I ) may break 
down as discussed in more detail in the Appendix of ref . |§ . 



5. Mixing Angle for HBT 

The BECF-s are frequently ||28|| but not exclusively |29|] parameterized by some 
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version of the Gaussian approximation. The out-longitudinal cross-term of BECF 
has also been discovered in this context recently [TB|]. In order to identify how this 
term may come about, let us assume that 

S c ,*{r,f],r x ,r y ) = H*(r)G*(r)) I*(r x ,r y ). (13) 

In Gaussian approximation one also assumes that 

F*(r) oc exp(-(r-7- s )7(2Ar, 2 )), (14) 
GM oc exp(-(r ] -r ]s ) 2 /(2Ar ] D), (15) 
h(r x ,r y ) oc exp(-( (r x - r X}S ) 2 + (r y - r yiS ) 2 )/(2Rl) ). (16) 

The corresponding BECF is given by a diagonal form as 

C(Ak; K) = l + K exp(-Q 2 At 2 - QffAtf - Q 2 Rl). (17) 

This diagonal form shall be transformed to an off-diagonal one if one introduces the 
kinematic relations between the variables Q T ,Q V and the variables QoutiQh- In the 
HBT coordinate system |27] one finds 

2R 2 out,LQoutQL),(l8) 

(19) 
(20) 
(21) 

-sinh[r /s ]) 2 Ar, 2 ,(22) 
(23) 

Note that the effective temporal duration, At* and the effective longitudinal size, 
T s Arj* appear in a mixed form in the BECF parameters SRl ut , R\ and R 2 out Ll and 
their mixing is controlled by the value of the parameter r] s . These results simplify a 
lotJJ in the LCMS, the Longitudinally Co- Moving System POfl, where /3l = 0: 

= PKcosh^Ar 2 + smh 2 [r ]s }r 2 Ar ] l), (24) 
R 2 L = cosh^^rfA^ + sinh^^lAr, 2 , (25) 
Rlut,L = -Asinh[r /s ]cosh[r /s ](Ar 2 + r 2 Ar ? 2 ). (26) 

Let us define the Longitudinal Saddle-Point System (LSPS) to be the frame where 
Vs( m t) = 0. In LSPS one finds that 

5R 2 out = f3 2 Arl (27) 
R 2 L = t 2 A V 1 + P 2 l AtI (28) 
Rlut,L = PtPL&r 2 . (29) 



C{Ak-K) 


= 1 + A* exp(-R 2 ide Q 2 side - R 2 out Q 2 out - R 2 L Ql - 


p2 

- rL side 


= R-l, 


p2 
n out 


= Rl + SRluv 


SRlut 


= /3 2 (cosh 2 [^]Ar 2 + sinh 2 [^]r 2 A77 2 ), 


Rl 


= ((3 L smh[i] s ] - cosh[?7 s ]) 2 r 2 A?7 2 + (/3 L cosh[?7 s ] - 




= (fa cosh[^](/5 L cosh [77,] - sinh[?7 s ]))Ar 2 + 




(P t sinh[r) s )((3 L smh[r] s ) - cosh[^]))r 2 A?7 2 . 
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Introducing Q = (3 t Q ut + PlQl and Q t = yQlut + Qlide ^ ne BECF can be rewritten 
in LSPS as 

C(Ak;K) = l + \*ex V (-AT>Ql-T*A v lQl-RlQ 2 t ). (30) 

Thus the out-long cross-term can be diagonalized in the LSPS frame |18| The cross 
term should be small in LCMS if r/f ,CMS << 1 i.e. if | y — Uq \« Ay ||. Since the 
size of the cross-term is controlled by the value of w s in any given frame, it follows that 
r] s is the cross-term generating hyperbolic mixing angle^ for cylindrically symmetric, 
longitudinally expanding finite systems which satisfy the factorization of eq. ([13]). 
The physical meaning of this hyperbolic mixing angle n s is illustrated on Figure 1. 



6. A New Class of Analytically Solvable Models 

For high energy heavy ion reactions, we model the emission function of the core 
with an emission function described in detail in ref.|J. This corresponds to a Boltz- 
mann approximation to the local momentum distribution of a longitudinally expand- 
ing finite system which expands into the transverse directions with a transverse flow, 
which is non-relativistic at the maximum of particle emission. The decrease of the 
temperature distribution T(x) in the transverse direction is controlled by parameter 
a, the strength of the transverse flow is controlled by parameter b. Parameter d con- 
trols the strength of the change of the local temperature during the course of particle 
emission||. If all these parameters vanish, a = b = d = 0, one recovers the case of 
longitudinally expanding finite systems with T(x) = Tq with no transverse flow, as 



discussed in ref . , if a = d = ^ b the model of ref . |T6[ is obtained. 

Specifically, we study the following model emission function for high energy heavy 
ion reactions: 



S(x;K)d 4 x = - o m t cosh(r] — y) H(j) x 
(2tt) 



/ K ■ u(x) M^) \ 
exp + — — dr T dr] dr x dr y . (31) 

V T ( x ) T ( X )J 

Here g is the degeneracy factor, the pre-factor m t cosh(r] — y) corresponds to the 
flux of the particles through a r = const hyper-surface according to the Cooper-Frye 
formula |31j and the four-velocity u(x) is 



/ 2 , 2\(l/2) / 2 , 2\(l/2) 

u(x) = Icosh^) l + 6 2 I^JM , b-, 6^, sinhfa) ll + b 2 x v « 



r o / r o T o \ T o 

+S),b\b r JL, 
2r J t t 



'cosh(77) (l + b 2 T I±ll) : b-,b^, smh( V ) (l + b 2 r l±J^L ) | . ,32) 




6 



which describes a scaling longitudinal flow field merged with a linear transverse flow 
profile. The transverse flow is assumed to be non-relativistic in the region where there 
is significant contribution to particle production. The local temperature distribution 
T(x) at the last interaction points is assumed to have the form 

1 + d 2 v - l° )2 ) , (33) 




T{x) T V 2t 2 ; V 2r, 



o 



and the local rest density distribution is controlled by the chemical potential fi{x) for 
which we have the ansatz 

MgO = tM) r\ + r 2 y (77 - y ) 2 

T(x) T 2R 2 G 2Ar] 2 ' 1 ' 

The parameters Rq and At] control the density distribution with finite geometrical 
sizes. The proper-time distribution of the last interaction points is assumed to have 
the following simple form: 

" (T) = ^ArV' 6 ^"^ 7 ' 2 ^- <35> 

The parameter At stands for the width of the freeze-out hyper-surface distribution, 
i.e. the emission is from a layer of hyper- surfaces which tends to an infinitely narrow 
hyper-surface in the At — > limit. 

The integrals of the emission function are evaluated using the saddle-point method 
p2 , |IT] , |17]]. The saddle-point coincides with the maximum of the emission function, 



parameterized by (T s ,rj s ,r XtS ,r y ^). These coordinate values solve simultaneously the 
equations 

dS _ dS _ dS _ dS _ q 
8t dr) dr x dr y 

These saddle-point equations are solved in the LCMS, the longitudinally Co-moving 
system, for v^ CMS << 1 and r XjS << r . The approximations are self-consistent if 
I Y — yo |<< 1 + Ar] 2 m t /T — Arf and (3 t « T 2 Ar^j (bR 2 ) which can be simplified for 
the considered model as fit = Vtl m t << {a 2 + b 2 ) jbj max(l, a, b). The transverse flow 
is non-relativistic at the saddle-point if (3 t « (a 2 + b 2 )/b 2 / max(l,a, b). We assume 
that At < tq so that the Fourier-integrals involving H(t) in the < r < oo domain 
can be extended to the — oo < r < oo domain. The radius parameters are evaluated 
here to the leading order in t x ^ s /tq. Thus terms of O(r XjS /T ) are neglected, however 
we keep all the higher-order correction terms arising from the non-vanishing value of 
7] s in the LCMS. 

For the model of eq. (|3TD the saddle point approximation for the integrals leads to 
an effective emission function which can be factorized similarly to eq. (|T3"|). Thus the 
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radius parameters of the model are expressible in terms of the homogeneity lengths 
Arj*,R*, At* and the position of the saddle point r] s i.e. the cross-term generating 
hyperbolic mixing angle. The saddle-point in LCMS is given by t s = tq, rj^ CMS = 
(y - Y)/{\ + A v 2 (l/A v 2 - 1)), r x , a = (3 t bR 2 J (t A v 2 ) and r y , s = 0. 



7. Geometrical vs. Thermal Length Scales for BECF-s 



The parameters of the correlation function are related by eqs. (|18|-|29D to the 
parameters R*, At* and T s Af]* which in turn are given by 



1 

J? 
1 

1 

At! 



1 1 

n G ^T 



A v 2 
1 



+ 



Arf T 



A 



11 + 



At 



cosh[^ s 
2 cosh[?7 s 



cosh [r} s 



(37) 
(38) 
(39) 



T 



Here the geometrical sizes are given by Rq, the transverse size, At], the width of 
space-time rapidity distribution and At, the duration around the mean emission 
time t s = tq. The hyperbolic mixing angle rj s ~ at mid-rapidity yo |J, where 
also the out-long cross-term [|16j vanishes. The thermal length-scales (subscript t) are 
given by 



Rn 



Tn 



a 1 + b 2 M t 



Arj. 



21 



At 2 



T, 



^0 

d 2 M t 



(40) 



K 2 



Kl 



The transverse mass of the pair is denoted by M t — y i\ 

These analytic expressions indicate that the BECF views only a part of the space- 
time volume of the expanding systems, which implies that even a complete measure- 
ment of the parameters of the BECF as a function of the mean momentum K may 
not be sufficient to determine uniquely the underlying phase-space distribution. 

It is timely to emphasize at this point that the parameters of the Bose-Einstein cor- 
relation function coincide with the (rapidity and transverse mass dependent) lengths 
of homogeneity^^ in the source, which can be identified with that region in coordi- 
nate space where particles with a given momentum are emitted from. The lengths 
of homogeneity for thermal models can be obtained from basically two type of scales 
referred to as 'thermal' and 'geometrical' scales. 

The thermal scales originate from the factor exp(— p ■ u(x)/T(x)), where u(x) is 
the four-velocity field. This is to be contrasted to the 'geometrical' scales, which 
originate from the exp(/j,(x) /T(x)) factor which controls the density distribution 0. 
Here p,(x) stands for the chemical potential. How do the changes of the temperature 
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in the transverse or temporal directions induce transverse mass dependent thermal 



radius or thermal duration parameters? See ref. [JXOj] for illustration. 

As a consequence of possible temporal changes of the local temperature, we find 
that the effective duration of the particle emission At* may become transverse mass 
dependent and for sufficiently large values of the transverse mass this parameter may 
become vanishingly small. The reason for this new effect is rather simple: Particles 
with a higher transverse mass are effectively emitted in a time interval when the local 
temperature (boosted by the transverse flow) is higher than the considered value 
for m t . If the local temperature changes during the course of particle emission, the 
effective emission time for high transverse mass particles shall be smaller than the 
effective emission time of particles with lower transverse mass values. 

For a more detailed analysis of the model the reader is referred to ref. || 3SJ, where 



it is pointed out that under certain conditions the parameters of the Bose-Einstein 
correlation function may obey an M r scaling: R S id e — Rout — Rl oc 1/VMt- 



8. Results for the Invariant Momentum Distributions 

The IMD plays a complementary role to the measured Bose-Einstein correlation 
function || |5], ||. Thus a simultaneous analysis of the Bose-Einstein correlation 
functions and the IMD may reveal information both on the temperature and flow 
profiles and on the geometrical sizes. 

For the considered model, eq. ([H|), the invariant momentum distribution can be 
calculated in such a manner, that the Cooper- Frye pre-factor m t cosh(^ — y) is kept 
exactly and the saddle-point approximation is applied to the remaining Boltzmann 
and proper-time factors, exp(— p ■ u/T(x) + fi(x)/T(x))H(r). This calculation yields: 

N iAp) = 7A3 (27rA^, 2 r 2 ) 1 / 2 (27ri?, 2 ) ^m t cosh(rr;) exp(+A^ 2 /2) x 

/AT 

exp(-p ■ u(x s )/T(x 3 ) + fi(x s )/T(x s )). (41) 
The quantities A^ 2 and % are defined as 

= 2 = T-2 + ITT COsh [^L Ws = 1 r^T 2JX 2 ( 42 ) 

Arj^ A?7 Z Ar/y 1 + A777 ^Vt 

and the modified saddle-point is located in LCMS at f s = t s = r , 777, f^rj = 
fttbR* I (tqAtit) and f^7 = 0. The modified radius and life-time parameters can be ob- 
tained by evaluating the R* and At* parameters at the space-time rapidity coordinate 
of the modified saddle-point, R~^ = R*(r)s CMS —> rfc) and Ar* = /\r*{r]l: CMS — > 77J). 
Thus the modified quantities (indicated by over-line) differ from the unmodified pa- 
rameters of the saddle-point approximation by the contributions of the Cooper-Frye 
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pre-factor. Note, that R* ~ i?* and At* ~ At* in the mid-rapidity region, where 
ris° MS ,rk « 1- 

In eq. fl¥T| ) the exact shape of the four-velocity field can be used, which is given 
in the first line of eq. fl32|). When evaluating /x(xj)/T(xj) in any frame, the invariant 
difference 1% — yQ CMS = rJ^ + y — yo should be used. 

The momentum distribution as given in eq. ([HP can be rewritten into a more 
explicit shape, which is more suitable for analytic study. This can be done if one 
neglects terms of Oij^rf /tq) in the exponent, which is the same order of accuracy 
which has been utilized for the solution of the saddle-point equations. Further, a 
term in the exponent (m t /To) coshf^j) is approximated by its second-order Taylor 
expansion, (m t /T )(l + 0.5r^ 2 ). These approximations yield 



9 , n -x— 2 2M/2 , n 



AMP) = J^3 (2ttA^ 2 t 2 ) 1/2 (2^) —l mt cosh(rj;) exp(+A^ 2 /2) x 
expWTo) -p(- 2( 1" 2 "+1|)) X 

where the geometrical contribution to the effective temperature is given by Toijnt) = 
T R 2 G /R^{m t ) = (a 2 + b 2 )m t R 2 g /tq and the fraction / is defined as / = b 2 /(a 2 + b 2 ), 
satisfying < / < 1. The following relations hold: 

Ay 2 (m t ) = A V 2 + kr? T {m t ) and ^ = — — — + (44) 

T» T + T G (m t = m) T 

Thus e.g. the width Ay 2 (m t ) is dominated by the longer of the geometrical and 
thermal length scale, in contrast to the HBT radius parameters which are dominated 
by the shorter of these scales. That is why the IMD measurements can be considered 
to be complementary to the BECF data. 



9. Analysis of Limitations 

The simple analytic formulas presented in the previous sections are obtained in 
a saddle-point approximation for the evaluation of the space-time integrals. This 
approximation is known to converge to the exact result in the limit the integrated 
function develops a sufficiently narrow peak, i. e. both 

Ar£(y,mt)«l and Ar 2 (y, m t )/r 2 « 1 (45) 

are required. These in turn give a lower limit in m t for the applicability of the formulas 
for the class of models presented in the previous section. These limits were studied 
analytically in ref. 0, here we analyze them numerically in the subsequent parts. 
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Figure 2: 

BECF parameters in LCMS, the IMD and the essential 'small' parameters. Solid line 
shows the results of numerical integration, dash-dotted line indicates the analytical 
results. Note that the transverse flow is switched off and T(x) = Tq as follows from 
a = b = d = 0. The longitudinal size is chosen to be small, At] = 0.3. On the 
last sub-plot, solid line indicates A77*, dash-dotted line stands for r XtS /r s and dashed 
line for rj^ CMS . These quantities are required to satisfy A77 2 << 1, r^ s /r s 2 << 1 

and r]^ CMS << 1. The other conditions, Ar^/r s 2 << 1 and Ai]1 << 1 are weaker 
conditions and thus not indicated on this and the subsequent plots. 



Compared to the condition (^), the condition of validity of the calculation of 
the invariant momentum distribution is less stringent, since one needs to satisfy only 
Ar]^ « 1. In the mid-rapidity region where NA44 data were taken, one has r) s ~ 
and for finite systems one finds m t » Tq (2 — 1 /At] 2 ). Note that this estimated 
lower limit in m t is extremely sensitive to the precise value of Arj in the region 
At] ~ 1/V2 « 0.7. For finite systems, the region of applicability of our results 
extends to lower values of m t than for infinite systems which were recently studied 



in great detail in ref. ||19|. An upper transverse momentum limit is obtained for the 



validity of the calculations from the requirement 

r x , s /T < 1 (46) 

or r 2 s /rQ << 1. This condition and the requirement br XtS /r Q < 1 has to be fulfilled 
simultaneously. Finally we note that the linearization of the saddle-point equations 
assumes that 

^ CMS < 1 (47) 

and for the IMD calculations a less stringent condition 777 < 1 was assumed. 

When comparing to data, detailed numerical studies may be necessary ]19[ to 
check the precision of the saddle-point integration. In the subsequent chapters, we 
present these tests for the analytic results given above and indicate the numerically 
found validity or violation of the conditions given in this Chapter. 



10. Numerical test for the a = b = d = case 

Let us start the numerical analysis with a set of parameters, that satisfy both 
conditions (^) and (|46|) . This is the reason for choosing a = b = d = and Arj = 0.3 
which guarantees that << 1 is satisfied. We start the analysis at mid-rapidity, 
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Figure 3: 

Same as Figure [| but with increased longitudinal size, A?7 = 1.25. In the \ow-m t 
region Ar)l > 1, which results in a slight increase of the numerically calculated R out 
component as compared to the analytical result. 



Figure 4: 

Transverse flow, transverse and temporal temperature gradients are switched on. 
Parameters correspond to a scaling limiting case at mid-rapidity. The numerical 
integration (solid line) is performed for the approximate transverse flow profile. Dash- 
dotted line indicates the analytical results. For the IMD sub-plot, dashed line stands 
for the simplified analytical result, eq. (f|3j), and dash-dotted line indicates the slightly 
more precise eq. (|4"T|). 



Figure 5: 

Same as Figure [|, but now the numerical integration was performed with the exact 
transverse flow profile. 



Figure 6: 

Same as Figure f|, but the longitudinal size of the system is increased. 
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and devote the last chapter to studies in the target and projectile fragmentation 
region. Note that the out-longitudinal cross-term |L6] vanishes at mid-rapidity. 

Due to the large number of model-parameters a rather comprehensive analysis is 
necessary and it is important to check the effects from various type of parameters. 

As one can see on Figure ^| the agreement between the analytical results and the 
numerically evaluated radius parameters and IMD is almost perfect. Note that for the 
analytical calculations of the HBT radii in LCMS eqs. (|l~8| - |29| , |37| - ^9l) were used, the 
IMD has been analytically evaluated both from the slightly more precise eq. (|43"D and 
from the more approximate eq. (41), both yielding results indistinguishable from the 
numerically integrated curve. For the numerical results on the HBT radii, we have 
evaluated the model-independent Gaussian radii in LCMS, as given in refs. [17], |l"6|| , 
for the model emission function of eq. (|3T|). For the numerical evaluation of the IMD, 
we have utilized eq. (§). 

How significant is the limitation given by eq. (f45|)? Can we apply the analytical 
results to systems created at 200 AGeV bombarding energy, where Arj ~ 1.5? Ac- 
cording to Figures |3] and ^| the analytical results provide a good approximation even 
in this case, (note that Figure ^|is evaluated for Arj = 1.5 for a,b,d ^ 0). The relative 
error on these figures is maximal around m t = 240 MeV for pions, being about 10 % 
for At] = 1.2, 20 % for At] = 1.5. These errors are characteristic for the R ou t and 
Riong radius parameters. The side radius component and the spectrum is obtained 
with much smaller errors. The relative error of the analytical results vanishes with 
increasing m t for this case. Thus the analytic approximations are reasonably good, 
sometimes excellent iia = b = d = and Arj < 1.5. 



11. Flow and Temperature Profile Effects 

The cases with non- vanishing parameters a, h and d are studied in this chapter. 
Figures |] and [| compare the analytical results with the numerical ones for Ar] = 0.7 
case, Figure |6| shows the same for At] = 1.5. The model parameters of these figures 
correspond to an approximate m t scaling for the radius parameters of the BECF. On 
Figure [5], the model was integrated numerically utilizing the exact flow profile, as 
given by the first line of eq. fl3"2"|). On Figures |] and the approximation given in 
the second line of eq. (^) has been performed before the numerical integration. In 
this case, the high m t limit of the analytical calculations corresponds to that of the 
numerical one, as it should because the radius parameters of saddle-point integration 
decrease with increasing m t and thus the precision of the saddle-point integration 
increases with increasing value of m t . 

In the region oim t — m < 200 MeV, the relative error of the analytical calculations 
for the HBT radius parameters is about 10 %. In case of the IMD the relative error 
refers to that of the slope parameter, which is approximately 10 % on Figures [| and 
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Figure 7: 

Numerical integration done with exact flow profile, case Rg < Rri^n) for pions. 



Figure 8: 

Same parameters as for Figure 0, but for kaons. 

It is worth comparing Figure f| with Figure [5[ The analytical curves on these 
figures are the same, but the flow profile eq. fl32|) is exact on Figure [|, while the 
approximation in eq. (^) was performed before the numerical integration on Figure 
|. The side and the longitudinal radius parameters are not changed significantly, 
however the high m% part of the spectrum and the out radius parameter is enhanced, 
increasing the deviation between numerical and analytical results to about 20 - 25 % 
in the kinematic region m t —m < 600 MeV. Note that the two analytic approximations 
for the spectrum, eqs. (f4l|) and (|43|) yield slightly different results, as indicated by 
the dashed and dash-dotted lines on the IMD part of the Figures, respectively. 



12. Pions vs. Kaons 

One expects that the precision of the calculation is increased for particles heavier 
than pions, because the thermal scales and thus the radius parameters decrease with 
increasing transverse mass. This is indeed the indicated by Figures |7| - |K]. 



When evaluating the numerical integrals for these figures, as well as Figures [LI] - [If], 
the exact flow profile has been utilized. 

The radius parameters are dominated by the geometrical ones on Figures [j] and 
resulting in a weaker m t dependence, while on Figures |9] and [H] the HBT radius 
parameters are dominated by the thermal scales since Rg > Rr(m) for pions, and 
the radius parameters satisfy approximate m t scaling. On these figures, the precision 
of the analytic results is nowhere worse than 20 %, in some cases especially for kaons 
and for the side and the longitudinal radius parameter at higher mt the relative errors 
become less than 5 %. 



13. Beyond the Limitations 

We have already seen that the conditions for the smallness of Arj* are automati- 
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Figure 9: Numerical integration done with exact flow profile, case Rq > Rri^n) for 
pions. 



Figure 10: 

Same parameters as for Figure E3, but for kaons. 



Figure 11: 

Results for kaons, using a parameter set where the limitation x s (m t )/ro << 1 is 
violated: x s (m t )/T > 0.6. 



Figure 12: Results for pions, using a parameter set where the limitation 
x s (mt)/T « 1 is violated. 



Figure 13: 

Another violation of the limitation x s (m t )/r << 1, using a = and 6 = 1. 



Figure 14: 

Numerical versus analytical results at y — 2. Note that mid-rapidity is located at 
yo = 3 and that the longitudinal size is given by Arj = 1.5. Note that At = fm/c. 



Figure 15: 



Same as Figure 0, but in the target fragmentation region y = 1 and for At = 2 fm/c. 
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Same as Figure [TR for y — 6. 



Figure 16: 



cally satisfied at higher m t thus the requirement (^) is not very important. On the 
previous graphs, we have shown figures where the other condition for the validity of 
the calculation, eq. ( [46] ) has been satisfied. In this chapter we violate this condition 
eq. ( f4"6] ) by a seemingly innocent change in the values of the parameters a, b, Rq 
and To. This way the saddle-point x s (m t ) deviates more and more from zero and 
it moves outside the region where the linearization of the saddle-point equations is 



allowed. Thus the analytical approximations break down as illustrated on Figures [11 
- [T~3|. The most serious defect is observable in the sub-plots indicating the IMD. The 
slope-parameter is very badly estimated by the analytical results if the condition (|46|) 
is violated. A good numerical interpretation of the condition (^) can thus be for- 
mulated as r XtS /r < 0.6. (Here we assume b < 1.) If r XtS /r < 0.6, the numerical and 
analytical results are in 10 -20 % agreement at mid-rapidity, however if t xs /tq > 0.5, 
the reliability of the analytic calculation is lost. The precision of model calculations 
for pions are very sensitive to whether this condition is well satisfied or violated. 

Note that our analysis has been performed with T = 140 MeV fixed. Lower values 
of To improve the agreement between the analytical results and the numerical ones, 
while larger values of To in general make the deviations larger in a given kinematic 
region. 

Surprisingly, the HBT parameters for kaons are well reproduced even if H^/tq 2 « 
0.5, except the R ou t radius parameter at high m t . For pions, the relative error on the 
radius parameters may reach 50 % for the side and the longitudinal momentum com- 
ponent, while it may reach 100% for the out radius component. Thus the condition 
of eq. (f46|) is substantial, it is mandatory to check it before concluding about the 
validity of the analytical approximations e.g. in a given data analysis. 



14. Departure from Mid-Rapidity 

Finally we turn our attention to LCMS radius parameters and IMD results off the 
mid-rapidity region. We use such parameters that the analytical approximations are 
precise to 10 — 20 % relative errors at mid-rapidity. For y ^ y , the out-longitudinal 
cross-term picks up non-vanishing values. This cross-term must vanish at m t = m, 
it develops a large modulus around m t — m = 50 MeV, which decreases fast with 



increasing value of m t , Figures [14] - [L6|. Surprisingly, the precision of the analytical 
approximations for side, longitudinal radius parameters and for the IMD is improved 
by moving one or two units away from the mid-rapidity. This is due to the fact 
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that r XjS , R*, Arj* and At* all decrease with increasing rj s which corresponds to 
increasing values of | y — y |. Figure [L^ indicates that the analytic calculations 
may remain reliable even if y = 1 for a source centered around yo = 3 and having a 
width A77 = 1.5, except the analytical result for the cross term in the low m f region. 
However, the analytical results for the cross-term are below the 20 % relative error 
limit if m t — m > 150 MeV or if At] < 0.7. With other words, the deviation of the 
analytical result for the cross-term from its numerically calculated value decreases 
very fast with increasing m t . Figures [14]- [16| indicate, that the analytical results for 
Rsidei Riong and the IMD can provide an approximation with as small as 5 - 10 % 
relative error, for R out with 15 % relative error even if the departure from mid-rapidity 
is substantial, | y — yo |~ 2Ar/. According to Figure the analytical approximations 
are reliable for m t — m > 150 MeV if y = 6, for lower values of m t the longitudinal 
saddle-point r]^ CMS is badly determined, due to the break-down of the linearized 
saddle-point equations. 



15. Summary and Conclusions 

We have presented a detailed numerical testing of the analytical approximations, 
as were given in || describing a three-dimensionally expanding, cylindrically symmet- 
ric, finite system. We find that the typical precision of the analytical approximations 
is characterized by a relative error of about 10 - 20 %, but in some cases much more 
precise agreement is found. Numerically, we have found that the analytic approxi- 
mations are reliable (within 10 -20 % relative errors) if the transverse position of the 
saddle-point satisfies r x s /r < 0.6 in some kinematic region. Numerically we find 
that this condition is the most sensitive among all the analytically found conditions 
for the validity of the analytical approximations, and a small violation of this condi- 
tion may result in large deviations between the analytical and numerical values. This 
corresponds to the break-down of the linearized saddle-point equations in this case. 

Surprisingly, we find that moderate deviation from mid-rapidity improves the 
agreement between the analytical results and the numerically found values for the 
HBT radius parameters and the invariant momentum distribution. This is due to 
the decrease of the transverse position of the saddle-point and the decrease of the 
effective source sizes with increasing deviation from mid-rapidity | y — y |. For 
large deviations from mid-rapidity, e.g. | y — yo \— 2Ar/, the analytical values for 
the cross-term and for the HBT radius parameters are found to be unreliable in the 
m t — m < 150 MeV interval. However, for larger values of m ( — m the HBT radius 
parameters, including even the cross-term, are reproduced correctly even at y = 6. 
This lower limit m t > 150 MeV corresponds to the requirement r]^ CMS < 0.9. 

In conclusion, we find that the analytical results are more reliable than expected 
before, since they are reliable in a fairly large rapidity interval. When comparing the 
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analytic results to data, one has to check if the transverse position and the space- 
time rapidity of the saddle-point is not too large in LCMS. The essential numerical 
conditions of consistency are found to be r x ^/r < 0.6 and r]g CMS < 0.9. 
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